Large-order Perturbation Theory for a Non-Hermitian "PT-symmetric Hamiltonian 
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A precise calculation of the ground-state energy of the complex PT-symmetric Hamiltonian 
H — p 2 + jx 2 + i As 3 , is performed using high-order Rayleigh-Schrodinger perturbation theory. The 
energy spectrum of this Hamiltonian has recently been shown to be real using numerical methods. 
The Rayleigh-Schrodinger perturbation series is Borel summable, and Pade summation provides 
excellent agreement with the real energy spectrum. Pade analysis provides strong numerical evidence 
that the once-subtracted ground-state energy considered as a function of A 2 is a Stieltjes function. 
The analyticity properties of this Stieltjes function lead to a dispersion relation that can be used 
to compute the imaginary part of the energy for the related real but unstable Hamiltonian H — 
p 2 + \x 2 - ex 3 . 

PACS numbers: 03.65-w,02.30.Lt,11.10.Jj 



It has been conjectured Q that the spectrum of the complex Hamiltonian 



H 



i\x s 



(1) 



is real and positive. Although there is no rigorous proof of this conjecture, it has been argued || that the reality 
and positivity of the spectrum is a consequence of the VT symmetry of H . (Recall that the parity operation acts as 
V : p — > —p and V : x — > —x and that the antiunitary time reversal operation acts as T : p — > —p, T : x — > x, and 
T : i — ► —i.) The notion that VT symmetry can replace the much more restrictive condition of Hermiticity has been 
studied in the context of quasi-exactly solvable quantum theories ||, new kinds of symmetry breaking in quantum 
field theory and complex periodic potentials ||. There have been many other instances of non-Hermitian VT- 
invariant Hamiltonians in physics. Energies of solitons in Toda theories with imaginary coupling have been found to 
be real m. Hamiltonians rendered non-Hermitian by an imaginary external field have been used to study population 
biology m and to study derealization transitions, such as vortex flux- line depinning in type-II superconductors ||. 

In this paper we study the large-order behavior of Rayleigh-Schrodinger perturbation theory for the ground-state 
energy of the complex 'PT-symmetric Hamiltonian (|l|). Note that this Hamiltonian describes a + 1 dimensional </> 3 
field theory, and recall that tfi 3 theories were the first quantum field theories in which the divergences of perturbation 
theory were studied |I(J . For the Hamiltonian (Q) we find that the perturbation series for the ground-state energy is 
divergent but Borel summable. Furthermore, by studying the numerical properties of the Pade approximants we infer 
that the (once-subtracted) ground-state energy considered as a function of A 2 is a Stieltjes function. This is a very 
strong result because it implies analyticity in the cut- A 2 plane and other properties. [It is surprising that this Stieltjes 
condition holds for a complex Hamiltonian such as (^); the proof that the once-subtracted ground-state energy of 
the conventional Xx 2N anharmonic oscillator is a Stieltjes function of A makes use of Hermiticity.] We then use these 
analyticity properties to establish a dispersion relation that yields the precise large-order behavior of the perturbation 
series. 

Let us consider the conventional Rayleigh-Schrodinger perturbation series about the ground state (E = \) of the 



harmonic oscillator H Q = p 



1^2 



The perturbed energy has an asymptotic series representation in powers of A 



because the perturbation x is an odd function of : 



25(A) 
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[We have chosen the form of Ho so that the perturbative expansion coefficients b n in (Q) are integers.] 

Using recursion formulas, we can easily generate as many terms as desired in this expansion. The coefficients b n 
alternate in sign, and their magnitude grows rapidly with n. The first 20 values are listed in Table 1. We have 
computed enough of the coefficients b n so that we can fit the leading large-n behavior as 



fi nn+l/2 
' (2^)3/2 



1-0 



(3) 



1 



Therefore, although divergent, the series in (||) is Borel summable [[Tljjl^]. Observe that if the factor of i were absent 
from the Hamiltonian (|l|) , then the perturbation coefficients b n would not alternate in sign and the perturbation series 
would not be Borel summable. 

We have performed a Pade analysis fll],[lj] on the divergent series for the once-subtracted ground-state energy 
[E(X) — ^]/A 2 . Using the first 46 perturbation coefficients b n , we find that for all real positive A 2 the diagonal 
Pade sequence P^(A 2 ) is monotone decreasing with increasing N, and the off-diagonal Pade sequence Pjv/ + i(A 2 ) is 
monotone increasing with increasing M: 

Pi <p}<b!<-..< pM+i < ■ ■ ■ < Pn < ■ ■ ■ < Pi < Pi < pS- (4) 

The results for A = 0.125 are shown in Table 2. If the inequalities in (Q) hold for all N and M and for all real positive 
A 2 , then it is rigorously true that [E(X) — \]/X 2 is a Stieltjes function of A 2 p2| . This means that [E(X) — \]/^ 2 
is analytic in the cut-A 2 plane, vanishes as |A 2 | — > oo, and is a Herglotz function of A 2 . [A function f(z) is said 
to be Herglotz if lmf(z) is positive (negative) when z is in the upper (lower) plane.] The fact that (Q) holds for 
< M, N < 23 provide strong numerical evidence that [E(X) — ^]/A 2 is a Stieltjes function. We stress that this is a 
much stronger result than merely saying that the divergent series @) is Borel summable. 

Furthermore, in addition to the inequality in (Q), the limits of the two Pade sequences appear to be identical. 
Therefore, we can extract values for the Pade summed energy from the two Pade sequences. The best estimate for 
the ground-state energy is obtained by averaging the last diagonal and off-diagonal Pade approximants. (To obtain 
an estimate of the ground-state energy from this average we multiply the average by A 2 and add |.) The results are 
shown in Table 3 for various values of the coupling A. Previous numerical calculations of the ground-state energy 
were obtained by direct numerical integration of the Schrodinger equation (see Ref. Q ) ; this technique gave a typical 
accuracy of about five decimal places. The agreement between the method of numerical integration and the Pade 
summation is excellent. Moreover, for A < the Pade technique provides an accuracy of more than ten decimal 
places. The agreement is better for smaller values of A, as is expected, because of a faster convergence rate of the 
Pade sequence. 



TABLE I. The first 20 perturbation coefficients b n 


in the expansion (g) of the ground-state energy for the complex 


"PT-symmetric Hamiltonian (|l]). 




n 


b n 


1 


11 


2 


-930 


3 


158836 


4 


-38501610 


5 


11777967516 


6 


-4300048271460 


7 


1815215203378344 


8 


-868277986898581530 


9 


464025598165231889260 


10 


-274145574452876905074540 


11 


177549419941607942489064216 


12 


-125174233315525265299874890500 


13 


95490636687662293430130201941400 


14 


-78410748996991270671939611723389320 


15 


68982408758305101330092396215438198608 


16 


-64750700102454900598854145411501140103290 


17 


64606224564767863138999679663986778514033420 


IS 


-68291871149169980983310351232642663615057109020 


19 


76244729314392095958565433992857306551429203990968 


20 


-89660576791390730762095201994590409692301843683859820 
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TABLE II 


. The diagonal and off-diagonal Pade sequences P$ (A 2 ) 


and Pjv+i(-^ 2 ) evaluated at A = 0.125. Observe the rapid 


convergence 


and note that the inequalities in (I 


M) are satisfied. 




7VT 
iV 




T-l JV 

Pn 


pJV 
" JV+l 







11.000000000 


4.739290085 


1 




7.039037169 


5.696806799 


2 




6.347866015 


5.947600655 


3 




6.168265727 


6.026389220 


4 




6.110857028 


6.054574069 


5 




6.089906566 


6.065678176 


6 




6.081499968 


6.070392205 


7 




6.077873385 


6.072516805 


8 
9 




6.076216002 
6.075421823 


6.073522627 
6.074018882 


10 




6.075025816 


6.074272525 


11 




6.074821510 


6.074406195 


12 




6.074712942 


6.074478558 


13 




6.074653729 


6.074518675 


14 




6.074620680 


6.074541394 


15 




6.074601848 


6.074554510 


16 




6.074590917 


6.074562214 


17 




6.074584462 


6.074566813 


18 




6.074580592 


6.074569597 


19 




6.074578237 


6.074571306 


20 




6.074576787 


6.074572368 


21 




6.074575882 


6.074573036 


22 




6.074575311 


6.074573460 



TABLE III. The ground-state energy for the Hamiltonian (jl|) for various values of the coupling A; the ground-state energy 
was computed by Pade summation and by direct numerical integration. The Pade sequences were computed for the once 
subtracted energy [E(\) — ^]/A 2 . The diagonal Pade energy refers to the energy extracted from the diagonal Pade sequence 
Pn (-^ 2 )i an< i the off-diagonal Pade energy refers to the energy extracted from the off-diagonal Pade sequence P$ +1 (\ 2 ). The 
best estimate for Pade energy is the average of the diagonal and off-diagonal values. 



A 


Diagonal Pade energy 


Off-diagonal Pade energy 


Pade energy 


Numerical energy 


0.015625 


0.50263 


0.50263 


0.50263 


0.50263 


0.03125 


0.50998 


0.50998 


0.50998 


0.50998 


0.0625 


0.53393 


0.53393 


0.53393 


0.53393 


0.125 


0.59492 


0.59492 


0.59492 


0.59492 


0.25 


0.71305 


0.71284 


0.71295 


0.71294 


0.5 


0.91445 


0.89035 


0.90240 


0.90026 


1.0 


1.40007 


1.05817 


1.22912 


1.16746 


2.0 


3.16075 


1.14032 


2.15053 


1.53078 
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The above Pade analysis provides strong evidence that the once-subtracted ground-state energy is analytic in the 
cut- A 2 plane. Thus, we can derive a dispersion relation in the expansion parameter A 2 to deduce the leading behavior 
of the imaginary part of the energy for negative A 2 . Physically, this means that we can compute the imaginary part 
of the energy (and hence the decay width) of the unstable ground state of the real Hamiltonian 

H =p 2 + \x 2 - ex 3 . (5) 

Note that the ambiguity in the choice of the sign of the coupling e corresponds to choosing the sign of i in (§). This 
has no effect on the decay width; the sign simply distinguishes the direction (left or right) in which the potential in 
(pi) is unstable. 

In the t = A 2 plane there is a cut along the negative t axis, and in the standard way |r^-|l5[] the b n coefficients are 
related to the discontinuity across the cut by the exact formula 

i r°° a D{-t) 



71- Jo * t H 

where D{—t) (t > 0) is the imaginary part of E(X) — |, evaluated with A 2 negative. From the growth estimate (||) we 
deduce that 

i_ 

D(-t)~-f_P=[l + 0(t)] (t-0+). (7) 

Thus, the leading contribution (for small e) to the imaginary part of the energy for the unstable ground state of the 
Hamiltonian (^) is 

ME(e)] ~ ^Z*E> (e -> 0+). (8) 

There are several ways to check this result. First, it agrees with a direct leading-order WKB calculation jlfj] of the 
imaginary part of the energy of the unstable ground state of the real Hamiltonian (g) . Second, applying the "bounce" 
method JlTj to the real unstable Hamiltonian (|j) we find that 

Im[£(e)] bouncc ~ c S^ 2 exp(-S ) (e -► 0+), (9) 
where the action So of the bounce solution is given by 



* =2 I**n*-**=h (10) 

and c is a constant (whose determination requires the computation of a fluctuation determinant). 

Finally, the answer in (^) is in agreement with the variational perturbation theory analysis in Ref. |lq| . In fact, 
Ref. jl8| contains a higher-order WKB expression for Im[i?(e)]. Inserting this higher-order WKB result into the 
dispersion relation (ph, we obtain a WKB-based prediction for the corrections to the leading-order growth of the b n 
coefficients given in (H): 



(tnn+l/2 / i 
I.WKB.,/ 1^+1 r L i _ 



169 44507 



120(71- \) 28800(n- \){n- |) 
9563539 189244716209 



1920000(n- h)(n- |) 8294400000(n - |)(n- §)(n- |) 



2/V'" 2 

42943442679817 



331776000000(n- \){n - |)(n - §)(n - |)(n - §) 

342541916236654541 
398131200000000(n- \){n - |)(n - f ){n - |)(n - |)(n- ^-) 

933142404651555165943 



(11) 



143327232000000000(n- \){n - |)(n - f )(n - |)(n - |)(n - ^)(n - f ) 
With these higher-order corrections, this growth estimate of the b n coefficients is spectacularly accurate. For example, 



4 



lWKB 

— = 1.00000000807. 

046 



(12) 



To conclude we note that the strategy employed here to relate the large-order Rayleigh-Schrodinger perturbation 
theory coefficients of a stable (and Borel-summable) problem to the imaginary part of the energy of an unstable 
(and Borel-nonsummable) problem is familiar from the quartic double- well potential H = p 2 + \x 2 + gx 4 , which is 
stable when g > and unstable when g < |is|JT^Jl5|| . The novelty in this paper is that we begin with a complex 
Hamiltonian H = p 2 + jx 2 + iXx 3 which, despite being non-Hermitian, nevertheless appears to be stable in the sense 
that it has a real and positive (and discrete) energy spectrum and a Borel-summable perturbation expansion for the 
ground-state energy. We can then relate the large-order perturbation coefficients to the imaginary part of the energy 
of an unstable state of the real but unstable Hamiltonian H = p 2 + jx 2 — ex 3 . It is interesting to note that the quartic 
case is relevant to the physics of instantons while the cubic case is relevant to 'bounces' in scalar field theories 

jl7j and to string perturbation theory |2l]] . 
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